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Recent signals for narrow hadrons containing heavy and light flavours are compared with quark 
model predictions for spectroscopy, strong decays, and radiative transitions. In particular, the 
production and identification of excited charmed and cs states are examined with emphasis on 
elucidating the nature of + and I + states. Roughly 200 strong decay amplitudes of D and D s 
states up to 3.3 GeV are presented. Applications include determining flavour content in rj mesons 
and the mixing angle in P and D wave states and probes of putative molecular states. We advocate 
searching for radially excited D* s states in B decays. 



I. INTRODUCTION 

The discovery of the L> s (2317) (0+) and L> s (2460) (1+) 
mesons 0, Q, with masses considerably lower than ex- 
pected in potential models^, stimulated a range of the- 
oretical activity. The current knowledge of charmed and 
cs mesons, compared to the expectations of Ref 3] is 
summarised in Figs. 1 1121 for the D and D s systems re- 
spectively. There is a significant amount of consistency 
between theory and experiment, interspersed with addi- 
tional states that do not fit well into such a classifica- 
tion, of which the D s (0 + ) and D s (l + ) are particularly 
sharp examples. Attempts to accommodate these states 
have invoked a variety of mechanisms. One interpreta- 
tion is that they are indeed qq 3 Pj levels, with their low 
masses being a realisation of chiral symmetry such that 
m(0+) - m(CT) = m(l+) - m(l~)_4]. An alternative 
is that they are multiquark or molecular configurations 
associated with the DK and D*K thresholds. One 
unresolved issue in the latter class of models is whether 
there are also further (0, 1) + broad cs states above D^*'K 
threshold, analogous to what appears to occur with light 
flavoured scalar mesons To help decide among com- 
peting interpretations, a coherent study of the dynamics 
of heavy-light hadrons is merited. 

A particular issue in testing these hypotheses will be 
to determine the 3 Pi — 1 P\ mixing angle for the ax- 
ial mesons. In the chiral symmetry picture^) this is 
implicitly assumed to be the ideal heavy quark limit 
(see Section Till B|) . In the molecular picture it is moot 
whether there is any simple mixing angle involving the 
D s (2460) and D s (2535) or whether a further axial with 
mass ~ 2.5GcV is called for. 

Subsequent to the above discoveries, SELEX0 re- 
ported a narrow state I? s (2632) seen in D s r] and DK 
with a branching ratio Z? s (2632) — > D s ri ~ QDK. The 
narrow width and the anomalous branching ratios (the 
two channels share the same quark flavours and phase 
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FIG. f: D Spectrum. The lines are predictions from Ref. [3j; 
the boxes are data 0. Both the Belle and Focus Do states 
are shown [2ll |. 



space favours the DK mode over the D s rf) led to sug- 
gestions that this state may be a tetraquark^ . Within 
the more conservative cs picture it was noted that the 
radially excited 2 3 Si is predicted to lie at ~ 2.73 GeV 
and that the presence of nodes in the wavefunction could 
lead to suppression of certain modes if the decay mo- 
mentum coincides with a node in momentum space 0- 
Such dynamics have been applied to the decays of ex- 
cited cc states ^3 an d also to light flavours with some 
success^!]; however, Ref. [3 found that such an explana- 
tion would only work if extreme values for the parameters 
were chosen and thereby concluded that the SELEX state 
might be an artefact. While we still agree with that con- 
clusion, it does raise the possibility that narrow states 
could in principle occur if their masses and decay kine- 
matics cause the momenta to coincide with nodes; this is 
one of the questions that we pursue in this survey. 

Independent of these tantalising possibilities, the tran- 
sitions among excited hadrons can determine their dy- 
namics and discriminate among models. With the pro- 
duction of charm in B decays, and the possible advent of 
charm factories at CLEO-c and GSI, it is timely to assess 
the landscape for heavy-light hadrons. 

Some highlights of the results are as follows: 
Significant production of D s (2 3 Si) and D s (3 3 Si) is 
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FIG. 2: D s Spectrum. See Fig. [UThe dashed lines are the 
DK and L»*Js: thresholds. 

predicted in B decays. Similarly, B decays may be a 
significant source of D\ mesons and can be used to test 
axial mixing. 

Decay ratios such as D*' — > D s r//D*' — ► Drj are useful 
probes of the flavour structure of the rj meson. 

The decays D*" — ► D s irj and D' sl r] test the putative 
molecular nature of the D s (2460). Similarly, the decays 
D'^ -> D s0 ry, D" -> A^if, and £>*" -> D sQ K* probe the 
structure of the enigmatic D s (2317). 

There may be considerable spectroscopic mixing be- 
tween the D s (2 3 Si) and _D s ( 3 Z?i) states which can be 
tested by measuring the transitions to DK and DK* 
from the vector D s states. 

El transitions such as 2 3 Si — > 1 Po,i are useful probes 
of the D\ and J5 s i mixing angles. 

Novel radiative transition selection rules are obtained 
in the heavy quark limit. 

Radiative transitions such as 1 + — > + 7, + — > 1~7j 
1 + — > 0^7 can test molecular models in the cs sector. 
For example, the width of the intermolecular transition 
1+ -> 0+7 ~ 17 kcV contrasts with the (9(1) keV rate 
predicted for cs states. 

Anomalous branching ratios of excited states, for ex- 
ample, D" s — > i?2^ being much larger than D*^ or 
DK*, may be used to probe the nodal structure of hadron 
wavefunctions. 

We proceed with a description of the strong decay 
model and the conventions concerning state mixing. This 
is followed by a discusson of the phenomenology of the 
strong decay results and their relationship to the heavy 
quark limit. The penultimate section concerns radiative 
transitions and highlights their utility as a diagnostic 
tool. All transition rates are contained in Appendices 
B and C; Tables llVl IVl and IVII for radiative transitions 
and Tables IVTlllXI VI for strong transitions. 



II. METHOD 

A. The Decay Model 

The 3 P model of strong decays assumes that qq pairs 
are created with vacuum quantum numbers|l3j|. Thus 
the interaction may be written as 

/ J 

where is a Dirac quark field of flavour /, nif is the 
constituent quark mass, and 7 is a dimensionless qqpair 
production strength. 

Recent variants of the 3 Pq model consider modifica- 
tions of the pair-production vertex [T^j , or assume that 
pair creation originates in a gluonic flux tube|l2l|. The 
latter is the "flux-tube decay model" , which in practice 
gives very similar predictions to the 3 Pq model. 

The model has been extensively applied to meson and 
baryon strong decays, with considerable success [n|, 0, 
[Tc| . The pair-production strength parameter 7 is fitted 
to strong decay data, and is roughly flavor-independent 
for decays involving production of uu, dd and ss pairs. 
A typical value obtained from computation of light me- 
son decays is 7 = 0.4[lll El 03, assuming simple har- 
monic oscillator (SHO) wavefunctions with a global scale, 
(3 = 0.35 - 0.4 GeV. The present work also assumes SHO 
wavefunctions but applies the formalism to a variety of 
heavy-light mesons; thus we have allowed the SHO (3 
values to vary according to the state (see Table [H] in Ap- 
pendix^). These values were obtained by equating the 
root mean square (RMS) radius of the SHO wavefunc- 
tion to that obtained in a simple nonrelativistic quark 
model with Coulomb+linear and smeared hyperfine in- 
teractions. Details are provided in the Appendix. 

In view of this it is appropriate to refit experimental 
data to obtain a new value for the coupling, 7. A total of 
32 experimentally well-determined decay rates have been 
fit with the model. Several variations of the decay model 
were also examined. Details and the final parameters and 
method are discussed in Appendix lAl 

B. Mixed States 

Heavy-light mesons are not charge conjugation eigen- 
states and so mixing can occur among states with the 
same J p that are forbidden for neutral states. These 
occur between states with J = L and S = 1 or 0. For 
example the J p = 1 + axial vector cn and cs mesons Di 
and D^ are coherent superpositions of quark model 3 Pi 
and 1 P\ states, 

|Di) = +cos(0)| 1 P 1 ) +sin(0)| 3 P 1 ) 

|Di) =-sin(</))| 1 P 1 )+co S (0)| 3 P 1 ). (2) 
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Quantifying the mixing pattern as a function of flavour 
will give information about the internal dynamics; little 
is known empirically at present. In the heavy quark limit 
Mq — > oo there is an explicit prediction for the mixing 
angle assuming that it is generated by spin-orbit inter- 
actions (see Section fill B|) . One of our aims will be to 
devise tests for determining this mixing in practice for 
heavy, but finite flavour masses. 

Mixing between S — 0, 1 states with the same J also 
occurs for 3 D 2 and 1 D 2 states: 

\D*) = + cos(0 D )\ 1 D 2 ) + sm((b D )\ 3 D 2 ) 
\D* 2 ') = -sinOMI^+cosOMI 3 ^) (3) 

Kokoski and GodfreyTlq find the angles (these su- 
persede those of Ref. PJ) <j)(lP cu ) = —26 degrees and 
4>(\P CS ) = —38 degrees, upon converting their mixing 
conventions to ours. Z?-wave mixing angles were not com- 
puted. 

Finally, the physical rj and rf are taken to be 

r) = + cos(9)—=(uu + dd) + s'm(9)ss 
v2 

rj = - s'm(9)-^= (uu + dd) + cos(9)ss. (4) 
v 2 

A typical quark model jj mixing angle is 9 — —45 de- 
grees. See Section fill El for further discussion on how to 
interpret the amplitudes for rj production. 

III. STRONG DECAYS 

We proceed with a discussion of OZI allowed strong 
decays of D and D s states. 

A. IS States 

D* decays are interesting because the open channels 
are so close to threshold that isospin symmetry breaking 
mass shifts become important. In particular D* + may 
decay to both D°ir + and D + n° , but the I? ± 7r =F mode of 
D*° is closed. The experimental D* + widths are in the 
ratio 



K* — » Kn is under-predicted by approximately 45% in 
amplitude. While this would be easy to correct by ad- 
justing the 3 Pq coupling, we choose to retain the fit value 
of Appendix IaI since it does well globally. We note that 
the decay model is tuned for SHO wavefunctions and de- 
cays of momenta of hundreds of MeV. Thus it is perhaps 
no surprise that the largest error seen in the 3 Pq model 
is seen in this extreme, near-threshold, decay. 

B. P- waves 

It is convenient to discuss heavy-light mesons in the 
jj coupling scheme. One has Lj = S\j 2 for the heavy 
quark which must combine with the light quark spin and 
angular momentum to form a total J p state. In the P- 
waves lj — P\j 2 and P^/2- Thus the j = 1/2 states form 
a doublet with J p = (0, 1) + while the j = 3/2 states 
form a J p = (1, 2)+ doublet. 

The relationship of these states to those in the LS 
coupling scheme can be determined once the heavy quark 
dynamics has been isolated and conventions have been 
fixed. We chose to employ the conventions of Ref. 1 171]. 
This reference also discusses the other conventions for 
the mixing angle that have appeared in the literature. In 
the heavy-quark limit a particular "magic" mixing angle 
follows from the quark mass dependence of the spin-orbit 
and tensor terms, which is 4>hq = —54.7° (35.3°) if the 
expectation of the heavy-quark spin-orbit interaction is 
positive (negative) • Since the former implies that the 
2 + state is greater in mass than the + state, and this 
agrees with experiment, we employ <f> = —54.7° in the 
following. This implies 

\Pi)hq = +^=\ 1 P 1 )- ] H\ 3 P 1 ) 

\P[)hq - +y| 1 P 1 ) + -^=| 3 Pi). (6) 

In practice the empirical mixing for the D and D s sys- 
tems is not yet known. Quantifying this is one of the 
challenges that we discuss here for both csand cqsystcms. 
For example, our decay model makes specific predictions 
for certain amplitude ratios: 



T(D*+ -+ D%+) A^P, - VPs) s = -^A( 3 P 1 -> VP S ) S (7) 

r(g*+ -> a**") = 2 - 21±0 - 57 - ( 5 ) v/2 

' ' = V2A( 3 P! -> VPs) D (8) 
The final state relative momenta are q — 0.0393 GeV i 



and q = 0.0379 GeV respectively, thus the form factor -4(2 J Si -► L PiPs) s = -^=A{2 3 S 1 -> 3 PiPs) s (9) 
is essentially unity, the ratio is dominated by the isospin 

factor and is predicted to be 2.28. A(2 3 S 1 ^ 1 P 1 Ps) D = V 7 2A{2 3 S l - 3 P 1 Ps) D (10) 

Absolute rates are given in Table Ivnl where one sees A( 3 D 1 — > 1 P 1 Ps) s — \/ r 2A( 3 D 1 — > 3 P 1 Ps) s (11) 
that the D* + widths are under-predicted by a factor of 

three. This is the largest error we have encountered which underpin the eventual extraction of the mixing an- 

with the 3 Pq model; for example the analogous decay gles. These relationships imply that the heavy quark Pi 
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state of Eq. couples to VPs in S- wave, whereas the P[ 
heavy quark state couples in Z?-wave. Thus one expects 
the D\ (Di ) to be broad (narrow) in the heavy quark 
limit. Similarly, Eq. 1111 implies that the the D\n (D[tt) 
mode will be large (small) in 2 3 S\ decays and the Ditt 
(D[ir) mode will be small (large) in 3 Di decays. 

Table IvTA give the predicted widths of the D\ and D[ 
states in terms of the mixing angle (c denotes cos <j)) . 
Equating these to the measured rates of 329 ± 84 MeV 
and 19 ± 5 MeV gives very good fits for mixing angles 
of 4> ~ —55° or 35°. The first angle is the solution for 
a broad D\ while the second is for a broad D[. Since 
these correspond to —54.7° and +35° respectively, good 
agreement is obtained with the heavy quark predictions, 
although distinguishing the two scenarios is impossible. 
This agreement may be tested by measuring decays with 
Dis in the final state. The decay tables indicate that 
the most promising such decays are D( 3 Di) — > Ditt and 
£>i7r, D*" -> Dxtt and D^tt, and D*" -► D X K and D[K. 

The situation for the D s i is less satisfactory because 
both the D s i and D' sl states are expected to be narrow 
due to the limited phase space available for the D* K 
channel. The heavy quark 3 Po model prediction for the 
D' gl width is 800 keV. Model uncertainties can be re- 
moved by measuring the ratio 



= 84 cos 2 +42 sin 2 6 S — 119 cos 6 S sin< 



r(£>^ -> D*K) 
T(D s2 -» D*K) 

(12) 

The D s i mixing angle may also be accessed through the 
decays D*" -► D sX K and D' sl K and D*" -> D sl r) and 
D' sl rj although only the D S \K mode presents a substan- 
tial branching fraction. 

The 1 + cs states can be directly produced by the W 
axial current in B — > DD s i or £>-Dsi[19|], thus B factories 
offer the possibility of studying these intriguing states. 
However, heavy quark symmetry suppresses the -P3/2 1 + 
decay constant^ so that production of the D' sl may 
be negligible. Furthermore there is no conserved vector 
current suppression of the scalar cs due to the differ- 
ent masses of the c and s in such transitions and hence 
the D s q may also be detectable. The relative production 
of these states in B decays will provide further insight 
into the relationship between the various cs states with 
J p = + ,l + . For example, the enigmatic D s q is pro- 
duced in D" and D*" decays to D sQ K and D s0 K* re- 
spectively; with the former having a branching fraction 
of approximately 17%. This is sufficiently large that it 
may worth looking for. One expects that these branch- 
ing fractions would be substantially lower if the D s o state 
were a molecule. 

In addition to the curious J p = + and 1 + states in 
the cs system, there are also questions in the charmed D 
states. A simple problem is the existence of incompatible 
candidates for the 0+ states at 2308 and 2407 MeVplj. 
The cd state may be produced directly in B decay as 
above, but with the additional penalty of Cabibbo sup- 
pression. It is also possible to produce the Dq in D" and 



D" decays. 

Finally, as noted in the Introduction, it has been sug- 
gested that the anomalously low masses of the Z) s (2317) 
and Z? s (2460) are consistent with breaking heavy quark 
and chiral symmetry 4]. This implies m(0 + ) — m(Q~) = 
to(1 + ) — m(X~ ). Presumably this relationship applies to 
S1/2 and pi/2 states and hence the mass of the broad 
Z?( s )i state should employed. For the D s system one ob- 
tains 349 MeV and 347 MeV for the left and right sides 
of the equation. Although not considered by the authors 
of Ref. 4] , in principle this relationship applies to the D 
system as well. The J — 1 mass difference is measured 
to be 347 MeV and the J = mass difference is 349 MeV 
if the lighter Belle Dq mass is used. Thus, if taken se- 
riously, this relationship supports the Belle results. We 
stress that it is important to test the applicability of this 
model through measurements of other observables, in- 
cluding the Pi mixing angle. 



C. 2S and ID waves 

All J p combinations of D* can be produced in decays 
such as B — ► D*iv 1 and the hadronic analogues, though 
transitions to states where the light degrees of freedom 
are in highly excited states will be suppressed by poor 
wavefunction overlaps and restricted phase space. The 
production of excited 1~ and 0~ should be feasible as 
they can be emitted from the W in weak transitions such 
as B — > D s D: the form factor that is relevant here is 
driven by the wavefunction at short distances. A simple 
quark model computation confirms that the wavefunction 
at the origin is not strongly suppressed as the principle 
quantum number increases, and the recoil momentum is 
sufficiently low that there is very little dependence on 
the B — > D transition form factor. Thus, the rate B — > 
D*(n)D depends primarily on the phase space and one 
finds that the relative emission rates are 



B -> D*D : B -> D* D : B 



D*"D 



1 



0.35 : 0.03 
(13) 

As the branching ratio b.r.(B — > D*D) and b.r.(B — > 
D*D*) are each approximately 1%0 one expects a to- 
tal production of 2 3 S"i (£>*') at roughly 1%. And even the 
doubly excited D* will have a branching fraction of order 
10~ 3 . The charmed analogue 2 3 S'i(D* / ) is Cabbibo sup- 
pressed and can be expected at the 10~ 3 level. As noted 
in the previous Section, the axial and scalar cs states 
can be produced in this way. For example, the 0~ D' s 
can also be produced; its decay is predicted to be almost 
entirely into D*K. We therefore advocate that these 
states should be sought at high statistics _B-factories. 

The signatures for the vector states are as follows. 
The cs 2 3 S X (£>*') decays dominantly to D*K(80%) and 
DK(15%) with traces of D s -q and D*T]. The D X K chan- 
nel is closed; to access the D\ this way requires the 3 3 Si 
initial state (see the discussion in the next section). The 



5 



charmed 2 3 Si(D*') decays dominantly to £>i7r(75%) and 
D*tt(10%) with traces of Dn and D*r\. 

The first excited vector D*' and D*' arise in either 
2 3 Si or l 3 Di configurations and in general there can be 
mixing between these. Such mixing will tend to shift the 
masses of the eigenstates away from the simple potential 
model values of Figs. 11121 The unperturbed masses for 
D* s of 2 3 S'i(2.73) or l 3 L>i(2.90) imply that one of these 
eigenstates will be kinematically forbidden to decay to 
D(IP)K while the other will be allowed. In the latter 
case there are interesting nodal effects, whose character 
will depend on the mixing angle. 

Similar remarks hold for the cq states and the decays 
to D(1P)it. However in this case the small pion mass 
implies that the D(\P)ir channel may be open for both 
initial states. The role of decays of these states in de- 
termining the 1 + mixing angles by decays to l + 0~ was 
discussed in Section lill Bl 

The + IP state is clearly Qpi/2, so there is no mix- 
ing problem with this final state in the transition from 
2S(0~) — > + 0~. Thus the + may be accessed via this 
transition from a D S (2S) produced in B decay. So given 
a D S (2S) emitted from the W current in B decay, one 
may access the + state by the above transition. How- 
ever, the analogous D' production is Cabibbo suppressed 
and it decays dominantly to D*ir; the D' s analogously 
decays to D*K. It would be interesting to study transi- 
tions to the analogous D s states and determine whether 
the -D s (2317) and D s (2460) are cs or other compounds. 
This would require the initial state to be 3S in order to 
be kinematically open. 

Mixing in the 3 D 2 — 1 P>2 system may be addressed in 
the heavy quark limit of the constituent quark model. 
Diagonalising the the spin-orbit interaction yields the re- 
sults 

M( 3 D 1 ) = M(D* 2 ) = M -~(H« o ) D (14) 
M( 3 D 3 ) = M(D;') = M + (H q so ) D , 

and a 3 D 2 — 1 D 2 mixing angle of </>£> = —50.76°. Thus, 
the heavy quark £)-wave states are 

\d;) hq = yfi 1 ^) - ^\ 3 d 2 ) 

\D* 2 ')hq = yfr^ + yfl 3 ^). (15) 

As with P-waves, the 3 Po strong decay model makes 
specific predictions for D-wave heavy quark decay am- 
plitudes which may be useful in interpreting the spec- 
troscopy. Some of these are: 

A{ 1 D 2 ^VPs) P = -^Ia( 3 D 2 ^VPs) p (16) 



A( 1 D 2 ~*VPs) F = +^A( 3 D 2 ^VPs) F (17) 
ACD 2 ^iPs) = (18) 

The second of these is an example of the 3 Pq selection rule 
forbidding such transitions among ggspin singlets[24|. 

We note that, in analogy with the P-waves, the am- 
plitude ratios above imply that the D 2 decays strongly 
in P wave while the D 2 decays only in F wave and is 
thus narrower than the D 2 . As with P-waves, this con- 
clusion agrees with spin conservation in the heavy quark 
limit. Unfortunately, the ability to distinguish the states 
is weakened by the many other decay modes which exist 
for these states. Nevertheless, the D 2 and D 2 have total 
widths which depend strongly on </>£> and measurement 
of any (or several) of the larger decay modes will provide 
(over) constrained tests of the model and measurements 
of the mixing angle. 

Finally we note that the transitions 1 D 2 — > VV and 
3 D 2 -> VV proceed in 3 P 2 , 5 P 2 , 3 F 2 , and 5 F 2 waves, but 
never share a wave. Thus there is no 1 D — 3 D mixing 
due to VV loops. 

The effects of wavefunction nodes can be seen in Figs. 
and Here the partial widths are plotted for fixed 
P as a function of the mass of the initial state. For the 
D*' , nodes would significantly affect the total width if 
the mass were roughly 3.1 MeV. Alternatively, nodes di- 
rectly affect the width of the D" by suppressing the D*K, 
DK*, and D*K* modes while enhancing the D 2 K mode. 
The figure indicates that a D" at 3.3 GeV would have a 
substantial branching fraction to D*K* while one at 3.4 
GeV would have no D 2 K mode. Clearly these effects 
must be accounted for in the phenomenology of heavy- 
light mesons. 
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D. 3S and 2D waves 

Similar remarks apply here as to the 25* waves with 
the bonus that phase space for decays into the IP states 
is open leading to potential measures of the axial mixing 



angle (Tables lXIl and lXiVfl . The challenge is to produce a 
significant sample of these 35 states in i3-decays. Our es- 
timate of the branching ratio b.r.(B — > D*" [D + D* + ]) ~ 
suggests that this may be feasible. Decays to D\X 
are predicted to be ~ 50%, which may provide a signifi- 
cant source of the axial charmed mesons and a measure 
of their mixing. Decays to D s in and D' ls rj may also be 
measured and used to test if the D s (2460) is the cs part- 
ner of the D s i(2535) or an independent state such as a 
D*K molecule. 

The decays of D" and D*" can access the enigmatic 
D s (2317) and Z3 S (2460). For example, comparison of the 
predicted rates for D*" — > D sl -q and D*" — > D' sl n will 
test the putative cs nature of the axial states. 

In the charm D system the decays D*" — > D\K, D[K 
and D2K measure the axial mixing angles. A robust 
prediction of the model is that the sum of the two axial 
decay modes should be roughly 7.8 times as large as the 
D2K mode. 

Finally, the status of £> (2308) and £> (2407) candi- 
dates may be tested by searching for them in the Dq-k 
and DqT] decay modes of the D" where they have a non- 
negligible branching fraction. 



natively both are permitted in 3S — > IS transitions. We 
note that the flavour flow in the D and D s systems probes 
the quark content of the 77's. Thus 77's are produced via 
their nn content in the cq system, whereas in cs decays 
it is the ss components that are involved. 

Thus a comparison of cq — ► n+cq and cq — > 
n+cq probes the nfi content weighted by phase space 
and form factor effects. A rather direct measure of the 
nn vs. ss content of the 77 can be obtained by comparing 

T (D*(2S) -» D{\S)n{nn)) : V (D*(2S) -> D s (lS)<n(ss)) 



or 



r (D(2S) -> D*(lS)n(nn)) : V (D S (2S) -> D*(lS)r)(ss)) . 

though the predicted branching ratios are small. 

The predictions in Appendix C for 77 and 77' production 
assume that these states are equally weighted mixtures 
of nn and ss . Hence with conventions for the n wave- 
function specified above, the amplitude ratio D s — > D s n 
and D — ► Dn is v^tan^. This is modified by different 
wavefunctions, different quark masses, and different me- 
son masses. Direct computation gives the factors shown 
in Table U these should be multiplied by \/2tan6> to ob- 
tain physical amplitude ratios. Determining several of 
these amplitude ratios will provide a constrained mea- 
sure of the mixing angle and the efficacy of the decay 
model. 



TABLE I: Amplitude ratios probing r\ decays with \/2tan( 
removed. 









D' s —>D*ri 


D"^Dii 


D*' -+D*7] 






1.78 


0.918 


1.09 


0.483 



IV. RADIATIVE TRANSITIONS 

Radiative transitions probe the internal charge struc- 
ture of hadrons and are therefore useful in determining 
hadronic structure. In particular they can help distin- 
guish possible exotic molecular or tetraquark state in- 
terpretations of the D s (0 + ) and D s (l + ) and determine 
mixing angles. 



E. Probing the 77 system 



A. El and Ml Transitions 



The emission of an 77 in 2S — > IS transitions is kine- 
matically allowed, while 77' emission is forbidden. Alter- 



El radiative partial widths are evaluated with the 
dipole formula 



7 



T^,i(nSLJ 



n'S'L'j' + 7) = - Cfi 6ss> 



m s Q + m 



r 



f) a\(nLJ\r\n'L'j') | 2 W 3 ^ 



(19) 



where Q and Q are the quark and antiquark charges in 
units of I e |, a is the fine-structure constant, lu is the final 
photon energy, Ef is the final state's total energy, Mi is 
the initial state's mass, and the angular matrix element 
Cfi is 



J p = 1 + D^K molecules has been assumed. Evaluat- 
ing this expression with SHO wavefunctions for the D^*\ 
K, and D s states and setting the scale of the D s molec- 
ular wavefunction using the weak binding relationship, 

gives 



Cfi = max(L, L')(2J' + 1) 



L' J'S 
J L 1 



(20) 



Wavefunctions were obtained from a simple nonrela- 
tivistic quark model which employs a Coulomb+linear 
central potential with an additional smeared hypcrfinc 
interaction. Tensor and spin-orbit terms are neglected. 
Results for El and Ml radiative transitions assuming qq 
structure are given in Appendix IbI 

Since to leading order El transitions are diagonal in 
spin, they select the 3 Pi component of the 1l and 1h in 
processes such as Di — > D*j. Thus measuring these rates 
yields a direct estimate of the P-wave mixing angle. The 
most promising processes are D[ ~ * P*7 w 800 cos 2 4» 
and D[ — + Z?7 ps llOOsin 2 ^, since both are large and 
involve the narrow D[. Prospects for using El transitions 
to measure <p s are less promising since the D s rates are 
all O(10) keV. 



and 



T(D s0 (mol) -> D*7) w 25 keV, (22) 



T(D sl (mol) ->■ D s j) re 30 keV, (23) 



r(D si (moZ) D*7) w 30 keV. (24) 



These predictions may be contrasted with the 0(1) 
keV results for the analogous El transitions of simple 
cs states reported in Table IV! 




FIG. 5: Molecular Radiative Transition. 



B. Molecular Probes 

The peculiar properties of the P/ S j(2317) raise the pos- 
sibility that this is either a cs state with substantial ad- 
mixture of the KD continuum pjj. a tetraquark state[2r|. 
or a DK molecule [3 • In the latter case it is suspected 
that the same dynamics give rise to a D*K resonance 
which may be identified with the D s (2460). 

The El radiative transitions 1 + — > 70~(cs) and + — > 
7l _ (cs) involve the overlap of molecular and cs wave- 
functions as shown in Fig. [SJ The amplitude is similar 
to one derived for the radiative decay of the X(3872)p7^ 
and is given by 



A 



1 1 2 

l v/23 



1 1 
713 



«0 



There is also the intriguing possibility of a radiative 
transition between the molecular D S (D*K) and D S (DK) 
states. In this case the rate is driven by a virtual D* — > 
.D7 transition, as shown in Fig. 



FIG. 6: Molecule-Molecule Radiative Transition. 



The resulting rate is given by 



d d pd s k(j) mo i(p)(j) D (k - U - p cu ^)(j)K(k + 7;- Psut: 



P, 



MolMol 



PcsT,) 



2 r 2' 
e*(9,A) 



ir(L>* 



► P+ 7 ;g >t )P + (^) + 
P°7;<Z*)Po(9*) (25) 



2^ 



where pij — (m^ — mj)/(mi + rrij), q is the momen- 
tum of the final state photon, and equal admixtures of 
the charged and neutral components of the J p = + or 



(21) 

where the width for D* — > D7 is evaluated at the photon 
momentum relevant to the 1 + — > r y0 + process, q t This 
result has been obtained assuming that the momentum 
space wavefunction for the molecular state is strongly 
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peaked at zero momentum, which will be the case for 
weakly bound states. The molecular form factors are 
given by 



F a (q) = J 



d 3 k 



3 ^: ) (fc)^: ) *(fc- 



m K 



m D + rriK 



q) (26) 



where the wavefunctions refer to the bound D*K and 
DK systems respectively and a is a channel index de- 
noting the D^+K or D < ~*}°K + components of the 
molecules. 

Weak binding implies that q* w mo* —mp « 140 MeV. 
This is reduced by the mass factor m^/(m£)±m^) = 0.2 
in the argument of the wavefunction in Eq. 1261 Thus the 
product is 28 MeV which is much smaller than the typical 
momenta in the weak binding potential, \/2^TdkEb ~ 
200 MeV. Thus the form factor may be neglected and the 
predicted rate is given by the average of the neutral and 
charged Ml transitions, D* — > Dj. Consulting Table 
IVII gives our final estimate T(l + (mol) — > + (moZ) 7 ) « 
17 keV, which may give a measurable branching ratio. 
This may be contrasted with the analogous cs transition 
predicted to be (0.26 cos <j> s + 1.4sin^> s ) 2 « 1.0 keV. 

Measurement of the total widths of the D s q and D s i 
will be required before the emerging data can be com- 
pared with these predictions 0. 



C. Radiative Transitions in the Heavy Quark Limit 



A(V -» 3 P 27 ) A =o 



> 



(32) 



Similar results hold for the decays involving the 1#. 

These results simplify in the heavy quark limit, reveal- 
ing some intriguing selection rules. In particular, in the 
heavy quark limit one has <fi = —54.7°, which implies 



- A(V -> 3 P o7 ) 
A(V — > 1h7)a=o 



A(V -> 1l 7 )a=o 
A(V -> 1l7)a=i 
A(V -> 3 P 27 )a=o 



\a(v 



A(V 



3 P 2 7)a=2 
3 P 2 7)a=i 



(33) 
(34) 



(35) 



The first of these corresponds to a selection rule which 
implies that the transition is pure El with M2 = (this 
is immediately clear because the process 3 Si — > 3 Po7 
has no M2 amplitude) . Hence the equality of amplitudes 
implies this is true for the 1l rate as well. These rela- 
tionships also imply that (apart from phase space correc- 
tions) 



r(2 3 Si - 1 L7 ) = 3T(2 3 Si -» 3 Po 7 ) 



(36) 



and 



As there are no charge conjugation constraints on the 
D and D s systems, radiative transitions from 3 Si can 
reach both 3 Pj and 1 Pi states. The most general trans- 
formation property of the transition at quark level is de- 
termined by noting that a positive helicity photon can 
change the projection of the spin or angular momentum 
of a quark by one unit. For radiative transitions between 
S and P levels the most general transformation for the 
current-quark operator is thus AL + + BS+ + CS Z L + , 
with unknown strengths A, B, C (these can be calculated 
in specific models but here we wish to make more general 
conclusions). The relative amplitudes for transitions to 
the various J p states are then driven by Clebsch Gor- 
dan coefficients. The procedure is defined in Refs.p^ 
and gives for the relative radiative amplitudes to or from 
(0,1,2)+ g]|: 



A(V 3 P o7 ) 

A(V - 1 l7 )a=o 

A{V ^ U 7 )a=i 
A(V -> 3 P 27 )a=2 
A(V -> 3 P 27 )a=i 



y/3 



A-C 



B 

71 



(27) 



-= ((A - C) sincj) + B cos $28) 
v2 

A _B = )^ + CcQs(j){2Q) 



s/2j V2 



A + C 

— (a + — 



(30) 
(31) 



r(ii 



3 5 l7 ) = r( 3 P ^ 3 5 l7 ). 



(37) 



Any deviation from these selection rules will measure de- 
viations from the heavy quark limit and the heavy quark 
mixing angle. 

The selection rules can be understood more imme- 
diately in the heavy quark jj basis. Only the light 
flavoured constituent can change its quantum state in 
the limit vtlq — > oo. In this case, the 3 Si may be rep- 
resented as Qs(i/2)3s(i/2) and the three independent P- 
states are combinations of Q s {x/2) and q in P1/2 and p 3 / 2 
with j z — ±1/2 or p 3 / 2 with j z = ±3/2. The three inde- 
pendent combinations of amplitudes in Eqs. I35l then cor- 
respond to the following: the transition si/ 2 — > P1/2 con- 
trols the identical amplitudes to the + and l + (si/ 2 Pi/ 2 ) 
states; the transition s 1 / 2 — > p 3 / 2 with j z — 1/2 deter- 
mines the second relation; and the transition s 1 / 2 — > p 3 / 2 
with j z = 3/2 controls the third relationship. 

With a large enough sample of 2 3 Si cs states, the ra- 
diative amplitudes to the cs 2 + and states can be 
compared with the selection rules to determine the mix- 
ing angle for the 1r state in the 3 Pi — 1 Pi basis. If the 
l + (2460) and + (2317) are the remaining states in the 
cs P-wave system, the same mixing angle should emerge 
when extracted from radiative transitions involving this 
pair of states. If O(10 3 ) events are required to measure 
radiative amplitudes, and each of these states is produced 
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with 6.r.(~ 10~ 3 ), then approximately 10 6 initial 2 3 Si 
mesons are required. Our estimates are that these arise 
at ~ 10~ 2 in B decays and so a suitable statistical sample 
should be accumulated at LHCb and other B-factories. 



V. CONCLUSIONS 

Excited D and D s states beyond the P-wave have not 
yet been identified. Moreover, mixing angles within the P 
and P-waves are not yet quantified. Determining these 
observables is an important task for spectroscopy and 
heavy quark physics, and forms a vital prerequisite for 
electroweak and CP violation studies. 

We expect that the emission of radially excited D s and 
D* will be significant in P-decays and can be anticipated 
at the 1% branching fraction level. Their decays give a 
potential source of P-wave cs states by which the 3 P— 1 P 
mixing may be measured. This in turn can determine the 
P3/2 — Pi/2 mixing pattern, which is important in testing 
models of the enigmatic P/ s (2317) and P s (2460) states. 
For example chiral models implicitly assume that these 
are pure pi/% configurations while there need be no simple 
mixing pattern in D^K molecular interpretations. The 
+ and l + (pi/2) states will also be emitted in P-decays 
via the W current. Their relative rates and the presence 
of one or two axial mesons in the data can also test the 
Pi/2 content and nature of the P s (2317) and P s (2460) 
states. Finally, radiative transitions also probe the quark 
structure of these hadrons and can assist in distinguishing 



molecular and cs assignments. 

Two body decay modes of excited D* states have nodes 
in simple nonrelativistic models. Establishing the reality 
of these is important as such nodes have been invoked 
to explain anomalies in the cc spectrum above charm 
threshold. The presence of nodes could in principle cause 
states to have narrow widths, though we find this unlikely 
in the decays considered here unless unfavoured values of 
parameters are employed. In practice the nodes for cer- 
tain decays tend to occur at energies where other chan- 
nels have opened, thereby restoring a canonical width for 
the state. However, the nodes are still implicitly mani- 
fested by virtue of the ensuing anomalous branching ra- 
tios; for example, nodes might suppress the DK chan- 
nel while allowing a £)(**) ifW decay, leading to relative 
rates that are an inversion of "phase space" expectations. 
Thus careful measurement of the relative branching ra- 
tios for 3 3 Si D s decays could prove to be a powerful tool 
for understanding gluodynamics in quark models. 
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APPENDIX A: DECAY COMPUTATION 
DETAILS 



1. Masses and SHO j3 Values 

The evaluation of the perturbative decay amplitudes 
requires mesonic wavefunctions. We follow tradition and 
employ SHO wavefunctions. Indeed, the model and ex- 
periment are sufficiently imprecise that computations 
with more realistic quark model reveal no systematic 
improvements p^j . The SHO wavefunction scale, denoted 
(3 in the following, is typically taken as a parameter of 
the model. However, since we seek to describe the de- 
cay of heavy quark states, it is preferable to fix the SHO 
scales to quark model wavefunctions. This was achieved 
by choosing (3 to reproduce the RMS radius of the quark 
model states. The resulting values are listed in Table HJ 

The meson masses used to determine phase space and 
final state momenta are listed below. 

Light meson masses: 

7T = 0.138, r) = 0.5477, p = 0.7758, w = 0.7826, K = 
0.495, K* = 0.8931, rf = 0.95778. 
D meson masses: 

D = 1.8694, D* = 2.0078, D 1 = 2.444, D[ = 2.422, 
D 2 = 2.459, D = 2.308 (Belle), D = 2.407 (Focus). 

Two experimental values for the scalar D meson mass 
are reported |2l| . Since these are incompatible we prefer 
to compute with both masses; leaving it to future exper- 



TABLE II: RMS Equivalent (3 Values (GeV). 



ft lj ,/ 


uu 


us 


ss 








1 Sio 


0.47 [0.4] 0.46 [0.4] 0.48 


0.43 


0.52 


0.71 




0.28 


0.32 


0.36 


0.37 


0.45 


0.66 


3 Pj 


0.26 


0.29 


0.32 


0.32 


0.37 


0.49 


0\Pi 


0.27 


0.29 


0.33 


0.33 


0.38 


0.50 


3 A/ 


0.25 


0.27 


0.25 


0.30 


0.34 


0.45 


O 1 ^ 


0.25 


0.27 


0.25 


0.30 


0.35 


0.45 




0.28 


0.29 


0.33 


0.31 


0.36 


0.48 


l 3 ^! 


0.24 


0.26 


0.30 


0.30 


0.35 


0.47 


2 1 S 


0.24 


0.25 


0.28 


0.28 


0.32 


0.41 


2 3 Si 


0.23 


0.24 


0.27 0.27 


0.31 


0.41 



imcnt to choose between the options. 
D s meson masses: 

D s = 1.9683, D* = 2.1121, D s0 = 2.317, D sl = 2.459, 
D' sl = 2.535, D s2 = 2.572. 
Theoretical masses were: 

D' = 2.58, D*' = 2.64, D" = 3.25, D*" = 3.31, 
D{ l D 2 ) = 2.83, D( 3 Di) = 2.82, D( 3 D 2 ) = 2.83, 
D{ 3 D 3 ) = 2.83, D' s = 2.67, D*' = 2.73, D" = 
3.24, Df = 3.29, D s ( y D 2 ) = 2.92, D S ( 3 D X ) = 2.90, 
D S ( 3 D 2 ) = 2.92, D S ( 3 D 3 ) = 2.92. 

These were obtained from the quark model used to 
determine the SHO scale or from Ref. 0. 

We set Di = Di(2444) and D[ = £>i(2422) since the 
latter is much narrower than the former. We also set 
D sl = L> s i(2459) and D' sl = D s i(2535) since both are 
narrow and the higher mass state is identified with the 
D' sl in the heavy quark limit. 

Finally, the quark model employed to determine the 
RMS (3 values and the radiative transition rates is a stan- 
dard color Coulomb plus linear scalar confinement inter- 
action with the addition of a Gaussian-smeared contact 
hyperfine term. The central potential is thus 

nr) = |c-^ + &r+|^5 CT ( r )S 1 -S 2 (Al) 
3 6 r 9mim 2 

where S a (r) = (cr/y/n) 3 e~ a r . The parameters were 
chosen to reproduce a broad range of open flavour masses 
and are C uc = -346 MeV, C sc = -319 MeV, b = 0.162 
GeV 2 , a = 0.594, and a = 897 MeV. Quark masses were 
taken to be m u = 0.33 GeV, m s = 0.55 GeV, and m c = 
1.6 GeV in both radiative and strong computations. 



2. Parameter Determination 

A variety of 3 Pq models exist. These typically dif- 
fer in the choice of weighting function used in the pair 
creation vertex, meson wavefunctions employed, and the 
phase space conventions. We shall restrict attention to 
the simplest vertex, which assumes a spatially uniform 
quark creation probability density. Possible phase space 
conventions include relativistic phase space (unit norm is 
used) : 
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TABLE III: 3 P Couplings. 



phase space 


Prms = 0.4 all = 0.4 


Rel. 
RPA 


0.485 ±0.15 0.417 ± 0.16 0.505 ±0.18 
0.214 ± 0.06 0.186 ± 0.07 0.228 ± 0.09 



(ps) = 2nk 



EbEc 



rriA 



(A2) 



where Eb is the energy of meson B in the final state. This 
can differ substantially from the nonrelativistic version: 



(ps) = 27rfc- 



m B m c 
m B + mc) 



(A3) 



especially when pions are in the final state. A third pos- 
sibility, called the 'mock meson' method, is employed by 
Kokoski and Isgur|12j: 



(ps) = 27rfc 



M B M C 
M A 



(A4) 



where Ma refers to the 'mock meson' mass of a state. 
This is defined to be the hyperfine-splitting averaged me- 
son mass. In practice, the numerical result is little dif- 
ferent from the relativistic phase space except for the 
case of the pion, where a mock mass of M r = 0.77 
GeV is used. The final possibility is referred to as 'RPA 
phase space' [2^ and postulates that the backward mov- 
ing Fock components of pseudo-Goldstone bosons (pions 
and kaons) contribute to decays. In the chiral limit the 
net effect of this is to multiply amplitudes containing a 
single pion or kaon by a factor of 2; if two Goldstonc 
bosons are present, the amplitude is multiplied by 3. 

We have investigated the efficacy of six models in de- 
scribing 32 well established experimental decay widths. 
These models all use SHO wave functions, either using 
a universal SHO width of = 400 MeV, the RMS- 
equivalent (3 values of Table |H1 or the RMS (3s with the 
exception of f3„ and (3k which are set to 400 MeV. The 
latter choice is an attempt to recognise that the lighter 
pseudoscalar states are Goldstone bosons and hence are 
likely to be larger than simple quark mode estimates. 
Relativistic and RPA phase space conventions have also 
been tested. We remark that the RPA and mock meson 
prescriptions yield similar results. 

The resulting best fit 3 Pq couplings and their errors are 
listed in Table EH As can be seen, the data and model 
are of sufficiently low quality that it is very difficult to 
distinguish the models (the models with (3^ = 0.4 obtain 
Sj/j « 38% while employing (3rms yields £7/7 ~ 29%). 



We henceforth adopt the relativistic phase space conven- 
tion with RMS (3 values determined as in Table [H] and 
set 7 = 0.485. 

The couplings required to reproduce experiment in the 
32 decay modes for the model used here are shown in 
Fig. [7| As can be seen, the large experimental errors 
preclude definitive conclusions. Nevertheless the model 
provides clear guidance over three orders of magnitude of 
predicted widths and over a broad range of quark flavours 
and meson quantum numbers. The specific decay modes 
are listed in Ref. 30] 
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FIG. 7: Equivalent Coupling vs. Mode. Po couplings re- 
quired to reproduce experiment for 32 decay modes. 



APPENDIX B: El AND Ml RADIATIVE 
TRANSITIONS 

Radiative transition rates based on Eq. ^| and the 
mixed states of Section Hi Bl are reported here. 



APPENDIX C: TABLES OF OPEN FLAVOUR 
STRONG DECAY MODES 

Strong decay rates and amplitudes are collected in Ta- 
bles IVHI to IXIVI In the following c and s refer to P- 
wave mixing in the D or D a sector c = cos (f> or cos <p s . 
Mixing angles in the D- waves are labelled c\ or s\\ thus 
ci = cos^D or cos 4>Ds- Estimates of decay rates con- 
taining mixing angles are given in terms of the theoreti- 
cal heavy quark predictions. Rates for modes containing 
the Belle and Focus Do are presented separately; when 
both may occur the differing possible total widths are 
separated by a slash. 
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TABLE IV: El Radiative Transitions in the D System 
mode q (MeV) F (keV) 
D+ Z; d* + j 408 51 
Dl -» D*°7 410 895 

(£) (J?) 

«) (£) 



'o T 7 



/311 
V223 



) (Ef) 



(£) («) 

-> D+7 175 9.4 

D"° -> D§7 175 163 

D+ -> D*+7 377 s 2 41 

L>? -> D*°7 379 s 2 715 

Di + -> D*+7 395 c 2 46 

Di° -> D*°7 398 c 2 819 

D*'+ -> D+7 209 s 2 9.5 

-> £)?7 209 s 2 164 

D*'+ -> Di'+7 189 c 2 7.0 

D*'° -> Di'°7 189 c 2 122 

D+ -> D+7 490 c 2 59 

D? -> D°7 493 c 2 1046 

-> £>+7 507 s 2 66 

Di'° -> £>°7 510 s 2 1154 

D'+ -> D+7 153 c 2 14 

D'° -> D?7 153 c 2 233 

D'+ -> D/+7 132 s 2 8.8 

£>' Di ,0 7 132 a 2 152 



TABLE V: El Radiative Transitions in the D s System 
mode q (MeV) F (keV) 



D a2 -> 73*7 


419 


8.8 


d s0 d* 7 


196 


1.0 


£>:' -» o s0 7 


382 


3.3 


d: # -> o S 27 


153 


1.2 


D s i -> D* sl 


323 


s 2 4.2 


D sl ' D*7 


388 


c 2 7.1 


D al D s7 


442 


c 2 7.3 


D s i' D s7 


504 


s 2 10.6 


d:' A,i7 


258 


s 2 3.2 


£>r - Dii7 


188 


c 2 1.3 




203 


c 2 5.4 


D' B Din 


132 


s 2 1.5 



TABLE VI: Ml Radiative Transitions 

mode q (MeV) r (keV) T (keV) PPG 

D *+ _^ jjj 1.5 ±0.5 

D*° -> D°7 137 32 < 800 

+ + /132\ / (0.7c + 1.8s) 2 
D 1 -^D ol ( 3? j ^ (011c + 28s j2 

132 \ / (3.2c + 2.0s) 2 

37 J V (0.47c + 0.31s) 2 

+ , + / (-0.56s + 1.4c) 2 

Ul ^ U " < \ 15 ) U-0.031S + 0.07c 



(-0.031s + 0.07c) 
, o nil\ / (-2.4s + 1.6c0 2 \ 
D 1 ^ A)7 { 15 ) {^Q.ns + O.OSc) 2 ) 
D* -> D s7 139 0.2 
D sl -> D s07 139 (0.26c + 1.4s) 2 

D' sl D s0 7 209 (-0.44s + 2.5c) 2 



TABLE VII: Open-flavor strong decays 



State Mode Y thy [T expt ] (MeV) Amps. (GeV~ 1/2 ) 



D* + D ,J 7T + 25 keV [64(15) kcV] 




1 j-, — 

Pi 


= -0.027 


D*+ D+7T 11 keV [29(6.8) keV] 




X A 


= -0.019 


D*° D'V 16 keV [< 2.1 MeV] 




J Pi 


= -0.021 


Do (2308) Dtt 316 [276(66)] 




'So 


= +0.635 


Do (2407) Dtt 283 [276(66)] 




% 


= +0.504 


Di D*tt HOc^ + 191s* - 230sc = 


272 [329(84)] 


3 Si 


= -0.338c + 0.478s, A Di = -0.153c - 0.108s 


D'j D*tt 191c 2 + 106s 2 + 240sc = 


22 [19(5)] 


"s 1 


= +0.504c + 0.357s, 3 Di = -0.099c + 0.141s 


D 2 Dtt 35 




L D 2 


= +0.164 


D*tt 20 




S D 2 


= -0.154 


Dj? 0.08 






= +0.013 


total: 55 [24(5)] 









TABLE VIII: Open-flavor strong decays, continued 



State Mode 


r thy \Yexpt 


] (MeV) 


Amps. (GeV- 1/2 ) 


D' D*7T 


27 




3 P = -0.145 


D (2308)tt 


41 




= +0.355 


D (2407)7t 


18 




= +0.416 


D*j? 


1 




3 P = -0.051 


total: 


69/46 






D* Dtt 


1 




1 P 1 = -0.025 


D*7T 


5 




3 Pi = -0.059 




0.4 




x Pi = +0.016 


D S K 


0.1 




x Pi = +0.010 


D*?7 


2 




3 Pi = -0.055 


DlTT 


11c 2 + 22s 2 


- 31sc = 33 


3 Si = -0.270c + 0.381s 








3 Di = -0.034c - 0.024s 


D' l7 r 


27c 2 + 14s 2 


+ 38sc = 1 


3 Si = +0.369c + 0.261s 








3 Di = -0.035c + 0.049s 


D 2 7T 


0.1 




5 Di = +0.031 


D S *K 


0.7 




3 Pi = -0.041 


total: 


44 
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TABLE IX: Open-flavor strong decays, continued 



State Mode r thy 




(MeV) 




Amps. (GeV _1/2 ) 




D( 3 Di) Dn 


73 








'Pi = -0.161 




D*7T 


45 








3 Pi = -0.141 




Dr? 


16 








1 Pi = -0.082 




D S K 


55 








x Pi = -0.164 




D*r? 


9 








3 Pi = -0.071 




Dl7T 


131c' 


! + 64s 2 


+ 183sc = 


0.2 


3 Si = +0.442c + 0.312s 
3 Di = +0.106c + 0.051s 




D'i7T 


61c 2 


+ 127s 2 


- 176sc = 


189 


3 Si = +0.288c - 0.409s 
3 Pi = +0.057c- 0.115s 




D27T 


7 








5 Pi = +0.108 




D S *K 


23 








3 Pi = -0.129 




Do 


74 








3 Pi = -0.207 






16 








3 Pi = -0.098 




D*o 


13 








1 Pi = -0.117, 3 Pi = 0, 5 Pi = 0.052, 5 Pi 


= 0.025 


D*w 


3 








x Pi = -0.063, 3 Pi = 0, 5 Pi = 0.028, 5 Pi 


= 0.011 


total: 


523 












Df d Di) 


53 








L F 3 = -0.136 




D*7T 


55 








3 P 3 = +0.155 




Dr? 


4 








1 F 3 = -0.041 




D S K 


4 








1 F 3 = -0.044 




D*r? 


3 








3 P 3 = +0.037 




DlTT 


3.5c 2 


+ 0.5s 2 


- 2.3sc = 


3 


3 P 3 = +0.071c- 0.026s 
3 G 3 = +0.013c + 0.009s 




D' l7 T 


0.6c 2 


+ 4.6s 2 


+ 2.8sc = 


2 


3 P 3 = -0.027c - 0.077s 
3 G 3 = +0.011c- 0.016s 




D27T 


6 








5 P 3 = -0.099, 5 G 3 = -0.010 




D S *K 


2 








3 P 3 = +0.034 




Dp 


15 








3 P 3 = +0.090 




Dw 


4 








3 P 3 = +0.050 




D*p 


99 








5 P 3 = -0.009, 1 F 3 = 0.343, 3 P 3 = 0, 5 H 3 


= 0.019 


D*cj 


27 








5 P 3 = -0.004, 1 F 3 = 0.188, 3 P 3 = 0, 5 H 3 


= 0.009 


Dr;' 


w 








1 F 3 = -5.5- 10~ 5 





total: 277 
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TABLE X: Open-flavor strong decays, continued 



State Mode 


r thy [r expt ] (MeV) 




Amps. (GeV _1/2 ) 


D 2 


D*7T 


93c? + 91s? + 9sici =87 




3 P 2 =+0.123ci -0.151si 
3 F 2 = +0.158ci +0.129si 




D (2308)tt 


0.8c? + 1.9s? + 2.4sici =0.4 




1 D 2 = -0.026ci - 0.041si 




D ( )(2407)7t 


0.9c? + 0.7s? + 1.6sici = 0.02 




1 D 2 = -0.034ci - 0.031si 




D* V 


13c? + 18s? - 22sici = 27 




3 P 2 = +0.077ci - 0.094si 
3 P 2 = +0.038ci + 0.032si 




DlTT 


0.7c? s 2 + 5.1s? c 2 + 7.6s?s 2 + 
3.8cisisc + 4.6cisis 2 + 12.5s?sc = 


0.1 


3 D 2 = -0.032cis - 0.087sic - 0.107sis 




D'jTT 


0.9c? c 2 + 9.8s? c 2 + 6.8s? s 2 + 
6.0cisic 2 — 5.0ciSisc — 16.3s?sc = 


8.4 


3 D 2 = -0.035cic + 0.096sis-0.115sic 




D 2 7T 


86c? + 115s? - 195sici = 197 




5 S 2 = -0.349ci + 0.427s! 
5 D 2 = -0.130ci +0.067si 
5 G 2 = -0.013ci -O.Ollsi 




D*K 


31c? + 45s? - 69sici = 73 




3 P 2 =+0.142ci -0.174si 
3 F 2 = +0.035ci + 0.029si 




Dp 

r 


75c? + 100s? - 122sici = 150 




3 P 2 = +0.183ci - 0.224si 
3 F 2 = +0.092ci + 0.075si 






25c? + 33s? - 41sici = 50 




3 P 2 = +0.106ci -0.130si 
3 P 2 = +0.051ci + 0.042si 




r 


50c? + 25s? = 33 




3 D 2 = -0.243ci, 5 D 2 = 0.172si 
3 P 2 = -0.025ci, 5 P 2 = 0.029si 






14c? + 7s? = 9 




3 Z> 2 = -0.133ci, 5 D 2 = 0.094si 
3 P 2 = -0.012ci, 5 P 2 = 0.014si 




D s0 K 


« 




1 D 2 = -0.007ci - 0.003si 




total: 


389c? + 437s? -409siCi + 8 








D*7T 


91c? + 93s? - 9sici = 96 




S D 2 = -0.151ci -0.123si 
3 P 2 =+0.129ci -0.158si 




Do(2308)tt 


1.9c? + 0.8s? - 2.4sici = 2.3 




1 D 2 = -0.041ci + 0.026si 




D (2407)7t 


0.7c? + 0.9s? - 1.6sici = 1.6 




1 D 2 = -0.031ci - 0.034si 






18c? + 13s? + 22sici = 4.5 




3 D 2 = -0.094ci - 0.077si 
3 P 2 = +0.031ci - 0.038si 




DlTT 


5.1c? c 2 + 7.6c 2 s 2 + 0.7s?s 2 + 
12.5c? cs — 3.8ciSics — 4.6ciSis 2 = 


1.2 


3 D 2 = +0.032sis - 0.087cic - 0.106cis 




D'jTT 


9.8c? c 2 + 6.8c? s 2 + 0.9s?c 2 - 

— 16c? sc — 6.0cisic 2 + 5.0cisisc = 


7.5 


3 D 2 = +0.035sic + 0.096cis - 0.115cic 




D 2 7T 


115c? + 86s? + 195sici = 3.9 




5 S 2 = +0.427ci + 0.349si 
5 D 2 = +0.067ci +0.131si 
5 G 2 = -O.Ollci +0.013si 




D S *K 


45c? + 31s? + 69sici = 3.3 




3 Z> 2 = -0.174ci -0.142si 
3 P 2 = 0.086ci - 0.035si 




Dp 


100c? + 75s? + 122sici = 26 




' i D 2 = -0.224ci - 0.183si 
3 P 2 = +0.075ci - 0.092si 






33c? + 25s? + 41sici = 7.9 




3 D 2 = -0.130ci -0.106si 
3 P 2 = +0.042ci - 0.051si 




DV 


25c? + 50s? = 42 




3 D 2 = +0.243si, 5 D 2 = 0.172ci 
3 P 2 = -0.025si, 5 P 2 = 0.029ci 




D*w 


7c? + 14s? = 12 




3 Z> 2 = +0.133si, 5 D 2 = 0.094ci 
3 P 2 = +0.012si, 5 P 2 = 0.014ci 




D s0 K 


w 




1 D 2 = -0.003ci + 0.007si 




total: 


437c? + 389s? + 409siCi + 9 
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TABLE XI: Open-flavor strong decays, continued 



State Mode 


r thy [r expt ] (MeV) 




Amps. (GeV" 1/2 ) 








D"(3.23) D*tt 


46 




3 P = +0.103 








Do(2308)tt 


72 




1 S = +0.154 








D (2407)7t 


63 




1 S = +0.156 








D*r; 


1.6 




3 P = +0.020 








D 2 7T 


36 




5 D = +0.124 








D S *K 


1.8 




3 P = +0.022 








Dp 


1.4 




3 P = -0.018 








Du 


0.6 




3 Po = -0.012 








D*p 


39 




3 P = +0.105 








D*LO 


13 




3 Po = +0.061 








D s0 K 


47 




1 S = +0.134 








D (2308)r? 


11 




1 S = +0.066 








D S K* 


2.2 




3 Po = -0.026 








D S *K* 


3 




3 Po = +0.034 










« 




5 Do = -0.004 








total: 


275/266 












D*"(3.31) Dtt 


53 




= -0.100 








D*7T 


59 




3 Pi = +0.112 








D77 


5 




1 P 1 = -0.033 








D a K 


14 




1 P 1 = -0.055 








D*r? 


4 




3 Pi = +0.030 








DlTT 


50c 2 + 62s 2 - 34sc = 


74 


3 Si = -0.095c + 0.135s 
3 £>i = -0.096c - 0.068s 








DiTT 


65c 2 + 53s 2 + 33sc = 


42 


3 Si = +0.134c + 0.095s 
3 Di = -0.070c + 0.099s 








D 2 7T 


35 




5 Di = +0.115 








D*K 


7 




3 Pi = +0.041 








Dp 


3 




3 Pi = +0.024 








Duj 


1 




3 Pi = +0.013 








D*p 


12 




1 P 1 = +0.012, 3 Pi = 0, 5 


Pi = 


-0.055, 5 Fi 


= 


D*a; 


5 




1 P 1 = +0.007, 3 Pi = 0, 5 


Pi = 


-0.034, 5 Pi 


= 


Dr)' 


0.5 




x Pi = -0.011 








D S K* 


0.7 




3 Pi = -0.013 








D s iK 


20c 2 + 34s 2 - 41sc = 


49 


3 Si = -0.085c + 0.120s 
3 Di = -0.039c - 0.028s 








DiiK 


27c 2 + 14s 2 + 36sc = 


2 


3 Si = +0.120c + 0.085s 
3 Di = -0.017c + 0.024s 








Dl»J 


4.3c 2 + 7.9s 2 - 9.8sc 


= 11 


3 Si = -0.041c + 0.058s 
3 Di = -0.016c -0.012s 










8.3c 2 +4.9s 2 +9.8sc 


= 1 


3 Si = +0.058c + 0.041s 
3 £>i = -0.014c + 0.020s 








D S *K* 


7 




x Pi = +0.011, 3 Pi = 0, 5 


Pi = 


-0.047, 5 Fi 


= 


D 2 j? 


0.6 




1 Pi = +0.017 








Do (2308) p 10 




3 Si = +0.066, 3 Di = 








D (2308)w 


3 




3 S"i = +0.038, 3 Di = 








D S „K* 


0.3 




3 S"i = +0.013, 3 Di = 








total: 


399/386 
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TABLE XII: Open-flavor strong decays, continued 

State Mode T thy [T expt ] (MeV) Amps. (GeV~ 1/2 ) 

D*K 261c* + 131s* - 369sc = 0.8 [< 2.3] 3 Si = 0.789e + 0.557s, J Di = -0.024c + 0.035s 

]%2 DK 27 L D 2 = +0.143 

D*K 3.1 3 D 2 = -0.069 





D B rj 0.2 


1 D 2 = -0.018 




total: 30 [15(5)] 






D*K 126 


3 P = -0.330 




F>*rj 0.5 


3 P = +0.043 




total: 127 




d;' 


DK 17 


x Pi = +0.090 




D*K 81 


3 Pi = -0.236 




D s ?7 2.6 


x Pi = -0.042 




DJjj 4.1 


3 Pi = +0.074 




total: 105 
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TABLE XIII: Open-flavor strong decays, continued 



State Mode 


Tthy [r expt ] (MeV) 


Amps. (GeV" 1/2 ) 




D s ( J Di) DK 


120 


"Pi = -0.205 




D*K 


74 


3 Pi = -0.181 






39 


1 P 1 = +0.129 






17 


3 Pi = +0.102 




DK* 


81 


3 Pi = -0.218 




total: 


331 






D s * 2 D*K 


155c? + 174s? - 92sici = 211 


S P 2 = +0.189ci - 0.232si 
3 F 2 = +0.173ci +0.141si 




D s *r, 


24c? + 35s? - 51sici = 55 


3 P 2 = -0.112ci +0.137si 
3 F 2 = -0.036ci - 0.029si 




DK* 


118c? + 165s? - 230sici = 258 


3 P 2 = +0.237ci - 0.290si 
3 P 2 = +0.089ci + 0.073si 




D (2308)K 0.04c? + 0.78s? + 0.34sici = 0.4 


1 D 2 = -0.006ci - 0.029si 




D s0 ?7 


0.07c? + 0.04s? + O.llsici w 


1 D 2 = +0.010ci +0.008si 




D*K* 


26c? + 13s? = 17 


3 P 2 = -0.207d, 5 P 2 = 0.147si 
3 P 2 = -0.008ci, 5 P 2 = 0.009si 




D (2407)K 


« 


1 D 2 = -0.003ci - 0.0055si 




D'jK 


w 


3 D 2 = -0.0009cic + 0.0026sis- 


- 0.0026sic 


total: 


323c? + 389s? - 373sici 






D* 2 ' D*K 


174c? + 155s? + 92s lC i = 118 


a P 2 = -0.232ci -0.189si 
3 F 2 = +0.141ci -0.173si 




D s *r, 


35c? + 24s? + 51sici = 4.0 


3 P 2 = +0.137ci +0.112si 
3 P 2 = -0.029ci + 0.036si 




DK* 


165c? + 118s? + 230sici = 25 


3 P 2 = -0.290ci - 0.237si 
3 P 2 = +0.073ci - 0.089si 




D (2308)K 0.78c? + 0.04s? - 0.34s lC i = 0.4 


1 D 2 = -0.028ci - 0.0006si 




D s0 ?7 


w 


1 P» 2 = +0.008ci - O.OlOsi 




D*K* 


13c? + 26s? = 17 


3 P 2 = -0.207si, 5 P 2 = 0.145ci 
3 P 2 = +0.008si, 5 P 2 = 0.009ci 




D (2407)K 


« 


1 P» 2 = -0.005ci + 0.003si 




d;k 


w 


3 D 2 = -0.0009sic + 0.003cis - 


0.003cic 


total: 


389c? + 323s? + 373sici 






D S ( 3 D 3 ) DK 


82 


L F 3 = -0.166 




D*K 


67 


3 P 3 = +0.168 




DsV 


4.5 


1 F 3 = +0.043 




D s *r, 


2.2 


3 P 3 = -0.035 




DK* 


14 


3 F 3 = +0.087 




D*K* 


52 


5 P 3 = -0.003, X F 3 = 0.294, 3 F 3 


= 0, 5 F 3 = 0.006, 5 H 3 = 


DiK 


« 


3 D 3 = -0.0004c - 0.0017s, 3 G 3 


w 


total: 


222 
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TABLE XIV: Open-flavor strong decays, continued 



State 


Mode 




| (MeV) 




Amps. (GeV 1/2 ) 


Di'(3.23) 


D*K 
DK* 


3.5 
1.9 
39 






3 P, = -0.030 
3 P = +0.024 
3 P = -0.100 




D (2308)K 47 






= +0.134 




D B 07J 


13 






= -0.073 




D*K* 


72 






3 P = +0.152 




Do(2407)K 35 






= +0.129 




D 2 K 


43 






5 D = -0.151 




total: 


219/207 








D s *"(3.29) DK 


9.6 






"Pi = -0.044 




D*K 


0.3 






3 Pi = +0.008 




D S 77 


0.4 






l Pi = +0.009 




D s > 


0.4 






3 Pi = +0.011 




DK* 


16 






3 Pi = -0.062 




D*K* 


117 






1 Pi = +0.040, 3 Pi = 0, 5 Pi = -0.180, 5 Pi = 




D s r/ 


0.5 






1 Pi = +0.012 




DiK 


18c 2 + 31s 2 


- 39sc = 


45 


3 Si = -0.082c + 0.117s 
3 Di = +0.035c + 0.025s 




D'jK 


33c 2 + 17s 2 


+ 44sc = 


1.6 


3 Si = +0.118c + 0.083s 
3 Di = +0.015c- 0.021s 




D 2 K 


6 






l Pi = -0.054 




D s i?7 


3.8c 2 + 6.9s 


2 - 8.9sc 


= 10 


3 Si = 0.041c- 0.058s 
3 Di = -0.015c- 0.010s 




~D'siV 


3.6c 2 + 2.4s 


2 + 3.5sc 


= 1.1 


3 Si = -0.045c - 0.032s 
3 Di = -0.015c + 0.022s 




total: 


208 









